Initial conditions and evolution of off-diagonal distributions 
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We briefly discuss the problem of specifying initial conditions for evolution of off-diagonal (skewed) parton 
distributions. We present numerical results to show that evolution rapidly washes out differences of input. 



1. Introduction 

Off-diagonal (or skewed) parton distributions 
provide important information about the nonper- 
turbartive structure of the nucleon [ |l[. They 
can in principle be measured in such processes 
as deeply virtual Compton scattering, diffractive 
vector meson production or diffractive high-py jet 
production. Just as for the ordinary (diagonal) 
parton distributions, the QCD evolution equa- 
tions [ ^] play an important role in the determi- 
nation of skewed parton distributions. As usual 
input distributions are required. However, this 
is more complicated than in the diagonal case 
since the skewed parton distributions depend on 
additional variable - the asymmetry parameter 
£ ~ p — p' . We use the symmetric formulation 
of Ji [ [l] in which the skewed parton distribu- 
tions are given by functions H(x, £) with support 
— 1 < ( x >£) < 1; see [ f° r more details. Here 
we discuss various ways to specify the initial con- 
ditions and show how the differences in the input 
distributions disappear on evolution. In partic- 
ular we illustrate the conclusion of [ |3) that the 
skewed distributions H(x, £), at small x and £, 
are fixed by the conventional diagonal partons. 

2. Constraints imposed on H(x, £) 

The distributions H(x,£,) have to fulfil several 
conditions with respect to the variable £. First, 
the time-reversal invariance and hermiticity im- 
pose the condition [ |l[ 

H(x,0 = H( X> -Q. (1) 



Thus H(x,!;) is an even function of £. The second 
condition states that in the limit £ = we recover 
the ordinary diagonal parton distributions 

H(x,0) = H diag (x). (2) 

The third condition is more complicated but has a 
simple origin, [[l], |j. The N th moment of H (x, £) 
is a polynomial in £ of the order N at most 

r l [N/2] 

/ dxx N - 1 H{x,i) = V A Nti e\ (3) 

J - 1 i=Q 

which also embodies condition (Q). Finally we im- 
pose continuity of H(x, £) at the border x = ±£ 
between two different physical regions, see [ |). 
This ensures that the amplitude of a physical pro- 
cess, described by skewed distributions, is finite. 
All these conditions have to be fulfilled in the con- 
struction of the input for evolution. They are, of 
course, conserved during the evolution. 

3. Initial conditions for evolution 

There are two equivalent ways to evolve 
H(x,^,/j, 2 ) up in the scale /s 2 . The first method, 
which is most appropriate at small x, £, uses the 
evolution of the Gegenbauer moments of H(x,£) 
and the Shuvaev transform [ 0, ^| to find the fi- 
nal answer in the a;— space [H|. In the second 
method the solutions are found after imposing 
initial conditions and numerically solving the evo- 
lution equations, directly in the x— space. Here 
we supplement the studies of [ |) by presenting 
results from the second approach. 
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Condition (|j|) is difficult to fulfil in order to 
specify initial distributions at a certain scale /Iq. 
It may be facilitated by the use of the double 
distribution !F{x, y) [ U 

H(x,£) = f dxdyT(x 1 y)S(x-(x + ^y)) 1 (4) 

where 1Z is the square \x\ + \y\ < 1 This pre- 
scription introduces a nontrivial mixing between 
x and £. Condition (Q) is guaranteed if F(x, y) = 
T(x,—y). We still have freedom to add to (Q) 
a function sign(£) D(x/£), antisymmetric in x/£, 
contained entirely in the ERBL-like region \x\ < £ 

0- 

The only problem left is to build in the ordinary 
diagonal distributions in the prescription . To 
do this we take 



T(x,y) = h(y)H d ^(£)/ 



l-\x\ 



dy'h{y>), (5) 



where different choices of h give, via (|J), different 
initial conditions for H(x,£). Three choices, 



m 

m = { a - 

sin(7ry 2 ) 



y 2 ) p+1 



p = 0,l, 



(6) 



are shown in Fig. 1, with H dla9 {x) given by 
GRV [ | at nl = 0.26 GeV 2 . The first choice 
gives simply the diagonal input H dla9 (x), inde- 
pendent of £. The second, with p = 0(1) for 
quarks (gluons), generates an input form simi- 
lar to that obtained from the model of [ |[ |(| in 
which the Gegenbauer moments of H(x,£) are £- 
independent. This property is conserved by the 
evolution. The exact form of the double distribu- 
tion in this case can be found in [ ||] . The last 
choice was selected so as to give an oscillatory 
input behaviour in the ERBL-like region. 

4. Discussion 

Fig. 1 shows the quark non-singlet, quark sin- 
glet and gluon skewed distributions for the three 
input models evolved up in /x 2 for £ = 0.03. The 
results for larger values of £ are qualitatively the 
same [[kJ. 

We see that the form of h(y) has the most in- 
fluence in the ERBL region. However even in 



this region evolution soon washes out the dif- 
ferences. Already by /i 2 = 100 GeV 2 the three 
curves are close to each other. In the DGLAP 
region, \x\ > £, the curves are almost identical 
while in the ERBL region, \x\ < £, they approach 
the asymptotic form for each particular skewed 
parton distribution. Such behaviour is especially 
important in view of the fact that one set of the 
curves is obtained from the pure diagonal input 
parton distributions. In this way we illustrate the 
result presented in [ |j| that to a good accuracy at 
small £, the skewed distributions il(cc, £; /z 2 ) are 
completely known in terms of conventional par- 
tons. Thus, to summarize, the nonperturbative 
information contained in the diagonal input par- 
ton distributions and particular features of the 
evolution equations for skewed parton distribu- 
tions are sufficient for their determination. 
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GRV ^=0.03 




Figure 1. The evolution to fi 2 = Q 2 = 4 and 100 GcV 2 of the quark non-singlet H NS , quark singlet H s 
and the gluon distribution H G starting from different inputs at /it 2 , = 0.26 GeV 2 for £ = 0.03 (indicated 
by the dotted vertical lines x = ±£). The dotted curves correspond to evolution from diagonal input, the 
continuous curves to the model of [3] and the dashed curves to oscillatory input in the ERBL-like region. 



